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Abstract 

We compare and discuss the respective efficiency of tfiree metliods (witli two variants for each of tliem), based 
respectively on Taylor (Maclaurin) series, Pade approximants and conformal mappings, for solving quasi-analytically 
a two-point boundary value problem of a nonlinear ordinary differential equation (ODE). Six configurations of 
ODE and boundary conditions are successively considered according to the increasing difficulties that they present. 
After having indicated that the Taylor series method almost always requires the recourse to analytical continuation 
procedures to be efficient, we use the complementarity of the two remaining methods (Pade and conformal mapping) 
to illustrate their respective advantages and limitations. We emphasize the importance of the existence of solutions 
with movable singularities for the efficiency of the methods, particularly for the so-called Pade-Hankel method. (We 
show that this latter method is equivalent to pushing a movable pole to infinity.) For each configuration, we determine 
the singularity distribution (in the complex plane of the independent variable) of the solution sought and show how 
this distribution controls the efficiency of the two methods. In general the method based on Pade approximants is 
easy to use and robust but may be awkward in some circumstances whereas the conformal mapping method is a very 
fine method which should be used when high accuracy is required. 

Key words: Two-point boundary value problem, Taylor scries method, Padc-Hankcl method, Conformal mapping, 
Polchinski's fixed point equation, Falkner-Skan's equation, Blasius' problem, Thomas-Fermi equation. 
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1. Introduction 



The use of expansions about the initial boundary (Taylor series) for solving an initial value (Cauchy) 
problem of an ordinary differential equation (ODE) is a well known procedure. Eventually an analytic 
continuation is needed to enlarge the limited range of validity of the series. It is less known that Taylor 
series may be used also to solve two-point boundary value problems (BVP) for ODEs. An old example is 

Email address: claude.bervillier@lmpt.univ-tours.fr (C. Bervillier''). 



Preprint submitted to Elsevier 



20 January 2013 



the Blasius problem [1] (a non linear ODE with specific conditions at the two boundaries of the domain 
P+ = [0, +00) of the independent variable, see section 3.5). In this case, the Taylor series is a Maclaurin series 
and, contrary to a Cauchy problem, one value is lacking at the origin (the unknown connection parameter ^ 
k) to make the series explicit. The Blasius method consists in first expressing the generic solution as a 
power series in the independent variable about the origin, the coefficients, c„ (fc), of which, depend on the 
connection parameter. In a second step, the missing value k* of k is determined by adjusting its value so 
that the sum of the series (or one of its derivative) matches the condition at the second boundary located 
at infinity^ . Of course, as in the Cauchy problem, the Blasius series has a finite radius of convergence and, 
at best, the "method can yield results of limited accuracy" (Weyl [2]). 

To improve the Blasius procedure, one must call for some analytic continuation of the series. But, unless 
one uses a computer and a symbolic calculation software, this is practically impossible because k is unknown. 
Moreover, with a nonlinear ODE, the possible appearance of spontaneous and/or movable singularities (the 
locations of which may depend on k) reduces, in an a priori unverifiable way, the magnitude of the radius 
of convergence of the Maclaurin series. This is why it is relatively recently that advances in the computing 
power have rendered practical the efficient use of several quasi-analytical methods for solving BVPs [3-9] 
(though the basic ideas are not always new). 

On general grounds, in a well posed BVP, one is seeking the unique solution that has some good analyticity 
properties, at least in The methods based on Taylor series for solving a BVP are implicit in that one 
attempts to determine the only value k* (the determination of an explicit solution in terms of the independent 
variable being left to a subsidiary step which often is obvious). For each method, there are two main variants 
for determining k* , according to whether the condition at the second boundary is explicitly imposed or not. 
We shall name them respectively "explicit" and "minimal" . 

"Explicit" procedure: In the first case, one constructs an explicit candidate solution (depending gener- 
ically on k) which is valid in the whole domain of definition I?+ of the solution sought. The value k* is 
then determined by imposing explicitly the condition at the second boundary (auxiliary condition). The 
explicit approximate solution is then obtained by substituting k* for k in the candidate solution. This 
procedure always works in principle (provided the candidate solution is well chosen) but in practice it 
may become very cumbersome due to the necessity of working with complete expressions. 

"Minimal" procedure: The second case is simpler since it amounts to imposing that the solution is 
merely analytically compatible with a particular class of functions having some (hopefully appropriate) 
analyticity properties in 2?+, such as generalized hypergeometric functions [3] or rational functions -the 
Pade-Hankel method [4-6]- or which satisfy a first-order polynomial ODE [7] or even the Taylor series of 
the generic solution itself -the simplistic method [10]. As in the "explicit" procedure, one first determines 
the specific function (hypergeometric or rational or etc... ) by the usual conventional matching with the 
coefficients of the Taylor series of the generic solution truncated at a given order M. Then, assuming 
that the function so constructed is analytically compatible with the solution of interest, one pushes the 
matching rules to the next order 1 all other things being equal. The auxiliary condition so obtained is 
used to determine k* . With such a procedure the requirement that the solution does satisfy the condition 
at the second boundary is not explicitly imposed. The "minimal" procedure is based on the hope that 
the solution sought has properties of analyticity similar to the function (or representation) considered. 
It may be efficient only when this solution is an isolated particular solution of the ODE defined in 2?+ 
(e.g. the envelope of a family of solutions having a movable singularity'*). In that case, imposing only 
specific conditions of analyticity in the domain may be sufficient to distinguish that solution among 
the others. Such occasions are not always realized, as, for example, with the Blasius problem which cannot 



^ We limit ourselves to problems where only one initial value is lacking (the connection parameter k) which is the value at the 
origin of the function, or of one of its derivative. It is to be determined in order to satisfy one condition at the second boundary 
located at infinity. 

^ In fact Blasius made use of a supplementary asymptotic expansion of the solution at the infinity point and matched the two 
expansions in a middle region of'D'^. This raises the supplementary question of the convergence of the asymptotic expansion. 
^ It may also belong to an infinite discrete set of solutions defined in D"'". 

* We explicitly show in section 2.3.1 that the Pade-Hankel method is equivalent to forcing the localization at infinity of a 
movable singularity (when it exists). 
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be solved with "minimal" proecdures, but "explicit" procedures work, see section 3.5 (this explains why, 
using the Padc-Hankcl method -the "minimal" variant of the Pad method- Amore and Fernandez have 
failed to determine the solution of that problem [11]). 

In this paper, after a brief presentation of the rough Taylor series method (and its "minimal" variant 
named the simplistic method Bervillier [10]), we discuss and compare the efficiency of two complementary 
methods of using Taylor series for solving a BVP. They rely upon the use of Pade approximants and of 
conformal mappings respectively. Of course, each of these methods has a "minimal" variant, named the 
Pade-Hankel method (Fernandez and coworkers [4-6]) in the first case and with no particular name in the 
second case [3,10]. Our choice is founded by their apparent easiness and efficiency relatively to other methods 
but also on their complementarity. As "explicit" procedures, the two methods are not new since they have 
already been proposed and used long ago to solve a BVP (e.g., see van Dyke [12, pp. 205-213] under the 
respective name of rational fractions and Euler transforms) . Only the smallness of the computing power of 
that time had limited their applications. More recent ("explicit") applications of the methods for solving a 
BVP may also be found in Boyd [13,14]. The use of "minimal" variants is less common. 

Our aim is twofold. Firstly we want to illustrate some conditions of application of the two methods 
with respect to the wide variety of domains of analyticity that one may encounter with a nonlinear ODE. 
The underlying idea is to try and propose some rules which may help in a priori determining whether the 
methods have a chance to succeed or not. Secondly we want to actually show that the two methods are 
complementary. In fact, the Pade method, when it works, is efficient, provided relatively low orders are 
concerned, when the desired accuracy requires the consideration of high orders of the series then the more 
sophisticated mapping method takes over. However, in turn, Pade approximants may be used to determine 
with a sufficient accuracy the two adjustable parameters involved in the mapping method. 

A short review on the use of Taylor series for solving a BVP for an example of a linear ODE (the 
eigenvalue problem of the anharmonic oscillator) may be found in [10]. Here we consider explicitly several 
nonlinear ODEs chosen for the different analyticity properties of their solutions in order to show (a small 
but significant part of) the variety of complexities that one may encounter. It is a commonplace to say that 
the disposition of the singularities in the complex plane controls the convergence of a Taylor series, it is 
less common to explicitly show this phenomenon in actual examples. The ODEs that we consider are the 
following (by order of complexity, each of them is accompanied by particular initial conditions not mentioned 
in this enumeration): 

- The Polchinski fixed point equation [15] in the local potential approximation (an exact, or nonperturbative, 
renormalization group equation, for introductory reviews see e.g. [16]). This is an easy example of second 
order ODEs having a single (most probably unique) non trivial solution defined in 2?"*" (it is the envelope of 
the general solution -a two-parameter family- having a movable singularity) . "Minimal" procedures (not 
the simplistic method however) succeed in determining this particular solution (due to its uniqueness). 
Although analytic in the (particular) solution of interest has singularities in the complex plane of 
the independent variable located exclusively on the negative real axis and the radius of convergence of its 
Maelaurin series is not small. The analytic continuations presently considered are then very efficient (see 
section 3.1). 

- The Thomas-Fermi equation [17,18] for the neutral atom, a second order ODE (and boundary conditions) 
with facilities similar to the previous case (general solution with a movable singularity, the solution of the 
BVP is the envelope of a family of such singular solutions). The main difference is that the singularities 
of the solution sought are no longer confined to the negative real axis but are arranged in the half-plane 
of the negative real axis (Ihs half-plane) . The convergence of the summation is then made more difficult 
than in the preceding ease, but the procedures considered (of any kind) work well again (see section 3.2). 

- The ODE which controls the flow of a third grade fluid in a porous half space [19,20]. With this equation the 
solution of interest has a Maelaurin series with a smaller radius of convergence and the general solution has 
an essential movable singularity at infinity (the leading part of the movable singularity is no longer a pole 
located at a movable point) whereas the solution sought (most probably unique) goes to zero. Procedures 
of any kind still work but they are not very efficient compared to the previous cases, it is preferable to 
consider the problem from the second boundary using a domain eompaetifieation: I?"*" = [0, -l-oo) — >■ [1, 0] 
via an essential change of definition of the independent variable (Ahmad [20]), for some values of the 
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parameters of the ODE, the new Maclaurin series then converges for k = k* (see section 3.3). 

- The Falkner-Skan flow equation. For the values of the parameters of this ODE that we consider, the 
singularities of the solution sought are disposed along an are of a circle around the origin (open toward 
the positive axis) with singularities that get in the half plane of the positive real axis (rhs half-plane) of the 
independent variable. This solution appears as the envelope of a two-parameter (particular) solution of the 
ODE with a movable singularity, it is (most probably) the unique solution defined in which satisfies 
the initial conditions. Procedures of any kind work, despite the slow convergence of the resummations (see 
section 3.4). 

- The Blasius problem. This is a particular case of the preceding ODE. The previous two-parameter (par- 
ticular) solution with a movable singularity disappears to leave room for a continuum of solutions defined 
in I?"*". The "minimal" procedures fail to determine a solution. Fortunately, the combination of a scaling 
property with the particular symmetry of the problem allows to transform the EVP into a simple Cauchy 
problem and makes the singularities of the solution sought disposed in the half plane of the negative real 
axis of the effective independent variable, consequently "explicit" procedures may be used with actual 
efficiency (see section 3.5). 

- The convex solution of the Blasius equation with more general initial conditions. It is the same kind of 
configuration as previously but the problem cannot be reduced to a Cauchy problem and the singularities 
of the solution of interest get in the right-hand-side half plane so that it is difficult to sum the Maclaurin 
series efficiently. Even the "explicit" procedures have difficulty to give an approximate solution (see section 
3.5.3). 

The organisation of the paper is as follows. We first shortly present the Taylor series method and the two 
analytic continuations that we consider. We then discuss each example of ODE independently. We finally 
conclude. In an appendix, we explicitly illustrate how to systematically implement a local expansion of a 
(general or particular) solution of an ODE. 



2. Taylor series method and analytic continuation 



2.1. Taylor series method 



This is an old and well known method for solving formally an ODE and several detailed presentations 
may be found (see for example [21]). We limit ourselves here to a brief introduction of its application to the 
EVP. 

Suppose that the function g (z) satisfies a second order ODE: 

F(5",g',5,z) = 0, (1) 

in which a symbol ' denotes a derivative with respect to the independent variable z. 

Because the general solution depends on two arbitrary constants, the BVP is to find (if it exists) the 
solution in a range ^ z G [a,b] which satisfies, e.g. the two following boundary conditions^ : 

(2) 

in which A and B are two given constants. It would be easier to determine the solution corresponding to 
two conditions given at only one boundary (the Cauchy problem), say, for example: 




gia) = A, 
g (a)^ k. 



(3) 



^ In fact, we shall have z S "D^ = [0, -|-oo) with the examples explicitly considered below. 

® The conditions could be, as well, expressed on the derivatives of g or mixed, provided that their number be equal to the 
order of the ODE. 
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In fact, the problem (2) amounts to solve the Cauchy problem (3) for an unknown (missing) value k* of 
k such that the condition g (b) = B be satisfied. If a unique solution exists, one may solve this problem 
numerically using, e.g., a shooting method and a Newton-Raphson procedure of adjustment. But presently 
one is rather interested in quasi-analytical methods based on the representation of the solution via Taylor 
series. 

To this end, one first Taylor expands the solution about the boundary a (Maclaurin series if a = 0): 

oo 

and determines the coefficients c„ in terms of A and k, such that the equation (1) be satisfied order by order 
in powers of [z — a). One gets that way the set c„ {A^ k) (with cq = A, ci = fc, etc.). It remains to sum the 
resulting series in order to satisfy the condition at the second boundary which determines the value k* of k: 

oo 

Y,c^{A,k*){b-aT =B. (4) 
Finally the global solution of the problem is (formally) given by: 

oo 

9* {z) ^^Cn {A, k*) {z - a)" . 

Of course the Taylor series method remains formal as long as one does not know whether the series 
converges or not. If it does, the method will only provide an approximate result since only series of a limited 
number of terms are considered in practice: 

M 
n=0 

At this stage we distinguish the "explicit" and the "minimal" variants of the Taylor method. 

2.1 A. "Explicit" procedure (Taylor method) 

The condition at the second boundary is explicitly imposed, namely: 

M 

Cn {A, k) {b-af^B, (6) 

n=0 

which provides the auxiliary condition to determine k*. This method is effectively used in section 3.3 where 
a useful change of function and a compactification of the initial EVP makes the resulting Taylor series 
convergent. But, in general, one encounters a difficulty when the second boundary goes to infinity (b — ^ 
+oo) because the Ihs of (6 ) goes to ±oo according to the sign of the last term. It is then better to switch 
to the "minimal" variant. 

2.1.2. "MinimaV procedure (simplistic method) 

The auxiliary condition is obtained by imposing that the next term of the truncated Taylor series vanishes: 

CM+i [A, A:) = . (7) 

That is equivalent to assuming that, for k = k* (ideally determined when M — > oo), the series converges 
within the boundaries and the hope is that the condition (7) will be sufficient to determine approximately 
k* for relatively small values of M . In the case where b ~ +cxj this would mean that the series of the desired 
solution has an infinite radius of convergence. In practice this is a so strong assumption that it is almost 
never true: the radius of convergence is most often limited by the presence of singularities in the complex 
plane of the independent variable z. Nevertheless it occurs sometimes that the simplistic method provides 
reasonable approximated estimations of k* (Margaritis et al [22]). For that two conditions are required: 
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(i) the solution sought is the unique analytic solution (or belongs to a finite set of separate analytic 
solutions). This is a condition required with any "minimal" procedure. 

(ii) the other solutions have movable singularities which must be (temporarily) located within the domain 
of convergence of the Taylor series whereas k is already close to k* ^ . 

Though the simplistic method does not work in general, it is so simple that it should be systematically 
tried before considering any other method. Unfortunately, it does not work with the set of ODEs considered 
below. 

When the Taylor series does not converge (almost always), one must call for analytic continuation methods 
such as Pade approximants or conformal mappings^ (see section 2.3). 



2.2. General remarks 



With the explicit examples of ODE and the methods considered in the present article, the auxiliary 
equations for k [similar to (6) or (7)] that we consider below are polynomial equations. The number of 
candidate values for k* (the zeros of the polynomial) equals the degree of the polynomial which grows with 
M . Even after having eliminating the complex zeros at a given order M, there is no possibility to distinguish 
the right value k* among all the remaining real zeros. But if the solution of the BVP is unique, one may 
expect to observe, when M grows, a clear convergence of one (and only one ^ ) series of (real) zeros toward 
the right value k* . When the method works, the convergent series may be distinguished alcady when M is 
small so that it is easy to locate and follow it for larger values of M (without having to determine all the 
other -uninteresting- zeros). Figure 4 of ref [10] illustrates how such convergent series of zeros may be easily 
located and followed. When no such series of zeros appears for small values of M, then the method works 
hardly or does not work at all. It is precisely one of the aims of this paper to illustrate how and why such 
situations may arise. 

The high degree of the polynomials to be treated makes it necessary to determine k* numerically so 
that the final expression of the solution is not completely put under an analytical form, the approximate 
method is called quasi-analytic. 



2.3. Analytic continuation methods 



2.3.1. Rational functions (Pade approximants) 

2.3.1.1. "Explicit" procedure (Pade method) Given the Taylor series (5) of finite order Af, one tries to 
approximate its sum by a rational function: 

, , P(z-a) 
Q(z - a) 

in which P (z) and Q{z) are polynomial functions of respective degree p and q. The coefficients of the 
monomials are unambiguously determined as functions of A and k from (5) so that: 

P{z-a)^Q{z-a) gu {z) + 0[{z- a)^'+^] , (9) 
Q (0) = 1 , (10) 
p + ci = M . (11) 



^ Actually, k* may be seen as the value of k for which the movable singularities are pushed to infinity, if those singularities are 
already outside the domain of convergence of the Taylor series when k is not close to fc*, the simplistic method cannot work. 
^ There are other possible analytic continuations. For example, the simplest one would be to consider the continuous analytic 
continuation method which consists in performing successive Taylor expansions each within its supposed circle of convergence 
until the second boundary be reached. But this method is not well adapted when the second boundary is located at infinity. 
^ This is still true if there is a discrete set of isolated solutions which satisfy the conditions, then one should observe several 
distinct convergent series of (real) zeros, e.g., see the eigenvalue problem of the anharmonic oscillator [10] where the generic 
eigenvalue parameter plays the role of the connection parameter. 

But with an unlimited accuracy if one uses a symbolic calculation software. 
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The representation (8) of gM (z) is also called the Padc approximant [p, q] of gM {z). It is a kind of analytic 
continuation of the truncated Taylor series (e.g., see Baker and/or Baker and Graves-Morris [23]). 

The solution of the boundary problem (1, 2) is then hopefully obtained by determining the values of k for 
which the condition at the second boundary is satisfied: 

P{b-a) _^ 

which is again a polynomial equation for k. 

The use of the explicit condition (12) at the second boundary to determine k* makes the procedure 
"explicit" . In principle g* (z) is unique and one should observe only one convergent series of zeros of the 
polynomial in fc as M grows. Notice that, for a given M , the pair (p, q) is not unique and may be chosen 
according to the expected behavior of g* (z) at the second boundary. In practice one may try several pairs 
{p, q) to obtain the best convergence of the series of zeros toward k* . 

The "explicit" procedure requires the explicit calculation (in terms of the unknown k) of all the monomials 
of the rational function (8) and this is generally very time consuming. Moreover, it may occur that the 
behavior of the solution at the second boundary is not compatible with a rational function and one must 
utilize some artifact to impose an efhcient condition [see, e.g. eq. (78)]. In addition, a Fade approximant 
may sometimes be ill defined [23] so that it cannot be an explicit representation of the solution and the 
method accidentally fails. 

The recourse to the "minimal" procedure allows to bypass these difficulties at least in the process of 
determining an approximate value of k* . 



2.3.1.2. "Minimal'' procedure (Pade-Hankel method) According to the presentation given in the introduc- 
tion, one imposes that the polynomials P and Q as defined at order M using (9-11) match again eq. (9) 
at next order (with M M + 1). The linear system of equations which then defines the coefficients of P 
and Q becomes overdetermined and leaves room for a definition of k via the vanishing of the determinant 
det(T) of a Toeplitz matrix Tij constructed with the coefficients c„ {A, k) [4-6]: 

Ti^j = Cq+^+i+,^j (i = 0,- • • ,(?; j = 0,- • • ,g) , (13) 

in which lu ^ p — q. 
Finally, the condition 

det(f) = 0, (14) 

provides a polynomial equation for k (the degree of which grows with M). 

Since the condition at the second boundary is not imposed explicitly, one is not assured (in particular if b 
is finite) that there is a unique solution and the set of zeros of (14) may or may not contain an approximate 
value of k* (see section 3). Let us show that, if the second boundary is located at infinity {b = -l-oo), the 
Fade-Hankel method has a particular significance. 



2.3.1.3. The Pade-Hankel method and the removal of a movable singularity Suppose that b = -I-cjo and 
that the generic solution of the ODE, corresponding to the initial conditions under study, has a singularity 
located at the movable point zq, then the solution g* (z), defined in is the envelope of that family of 
(singular at zq) solutions. This envelope may be seen as the limit zq — > +oo of the family. Consequently, to 
determine fc*, one may proceed as follows. 

Instead of imposing that the ratio P/Q should also reproduce the series at next order gM+i (z) one requires 
that the rational approximation has a singularity (a pole) located at some zq, i.e.: 

Q(^o) = o, 

this induces a modification of the matrix (13) which then reads: 

^^The matrix Tij may be rewritten under the form of a Hankel matrix Hi i by a redefinition of the index I = q — j . The 
reference to a Toephtz matrix is more convenient with regards to the discussion presented in section 2.3.1.3. 
Depending on the unknown connection parameter k. 
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Ti,j =c,+^+i+j_j (i = 0, • • • ,g - l;j = 0, • • • ,g) , 

By imposing that det(r) = one obtains a polynomial equation for k which depends on zq. When zo 
becomes very large, det(T) is dominated by the largest power of zq (zq) so that, when zq — >■ +oo, the 
effective auxiliary condition asymptotically reduces to the vanishing of a minor of the matrix T which is the 
determinant of the following reduced matrix: 

ft,j = Cg+^+i+.^-j (i = 0, • • • , g - 1; j = 0, • • • ,q-l) . 

That is nothing but the Tocplitz matrix (13) considered at the preceding order (g — > g — 1) but with 
w — > w + 1. 

Hence the Pade-Hankel method (for determining k*) is closely related to the process of sending a movable 
pole at infinity. Its efficiency thus depends on the existence of such a movable singularity. 



2.3.2. Conformal mappings 

In the following we assume that the two boundaries are a = and b = +oo. Let us consider a continuation 
of the independent variable z in the complex. 

The conformal mapping method for solving a BVP has been introduced by Bervillier et al [24] to which 
we refer the reader for more details. 

If one knows that the solution (i.e. when k — k*) is analytic in the right-hand side interior of an angular 
sector as drawn in figure 1 of [24], then one can effectuate the following conformal transformation (R > 0, 
a > 0): 

(l + z/i?)^/"-l 

z^w = — , (15) 

(1 + z/i?)^/" + 1 

which maps the region of the z-plane formed by the interior of the angular sector the vertex of which is 
located on the negative real axis at — i? and of angle cm 12 measured in radians into the interior of the 
unit circle centered at the origin of the ui-planc so that z ~ +cxd corresponds to w = +1 (whereas z — 
corresponds to w = 0). 

The transformation (15) is better understood as the result of the following sequence of tranformations: 

(i) z — Zi = z/i?, 

(ii) zi — >■ Z2 = zi + 1, shift of the vertex of the angular sector onto the origin, 

(iii) Z2 — >■ Z3 = (z2)^^", mapping of the interior of the angular sector onto the positive half plane, 

(iv) Z3 — >■ u; = (z3 — l)/(z3 + 1), mapping of the positive half plane onto the unit disc \w\ < 1. 

This conformal mapping is merely a generalization of the Euler transform which corresponds (sec, e.g. 
van Dyke [12, p. 208]) to set a = 1 and i? = 1/2 (the rhs half plane with boundary located at z = —1/2). 

If the interior of the angular sector is a region of analyticity of the original function g (z) then the series 
(5) converges there. Under the analytic continuation this convergence is conveyed to the whole unit disc for 
the series: 

M 

9M M = Ui {A, k) w' , (16) 

i=0 

obtained by expanding, within the original series (5), the relation inverse of (15): 

'1 + ' 



R 



1 



1 



(17) 



Before looking at the way one may determine the values of R and a, let us define the two variants of the 
mapping method. 



2.3.2.1. "Explicit" procedure The condition at the second (here infinite) boundary may be imposed by 
setting w = 1 in (16): 
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AI 

.9A/(1) (18) 

4=0 

2.3.2.2. "Minimal" procedure The condition (18) indicates that the mapping method is naturahy an "ex- 
phcit" procedure. It may be used also formaUy by simply imposing that the last term of the series (16) 
vanishes: 

UM {A, k)=0. (19) 

As in the case of the simplistic method (see eq. [7]), one forces the convergence of the series, but this 
time the series of the solution of interest actually converges provided the conformal mapping has been 
well chosen. This gives a simple and light version of the mapping method for determining k* however the 
gain compared to the use of the "explicit" procedure is not as large as that one may observe with Pade 
approximants. Moreover, as mentioned previously in section 2.3.1, the "minimal" procedure may not work 
at all (see section 3). 

2.3.2.3. Choice of the parameters R and a The determination of the parameters a and R depends on the 
distribution of the singularities of g* (z) in the complex plane which is a priori unknown. In practice one 
may try several values of a and R beginning with sufficiently small values so as to get a convergent series 
of zeros for k in (18) or (19). Then this rough estimate of k* may be used to determine a better value of R 
by the d'Alembert or the Cauchy rule [3,10] because, most often ^"^ the maximal value of R is the radius of 
convergence of the series (5) for k = k*. 

Actually the radius of convergence of the Taylor scries is fixed by the position in the complex plane of the 
independent variable of the singularity the closest to the origin (remind that a = 0). If one assumes that 
this position is known to be — x + iy, then one has: R = y^x"^ + y'^ and a ~ 2 arcsin (|y|//i) /tt (in which 
h= .yiR + x)^ +y^). 

We shall show that Pade approximants may be used to determine the values of R and a (in a way similar 
to that of Andersen and Geer [25]). 

3. Explicit examples of ODEs 

In the following we shall generically call g* (z) the solution of the BVP of interest, it corresponds to a 
(hopefully unique) value fc = fc* of the connection parameter which will be the "missing" value of g or of one 
of its derivative at the origin. The assumed domain of definition of g* (z) is P+ = [0, +oo), but sometimes 
solution defined in = [0, — oo) will be also encountered. 

3.1. The Polchinski fixed point equation in the local potential approximation 

The problem is to find the potential V (x) — V* (x) which satisfies the following second order ODE and 
boundary conditions: 

V" - v'^ - ^xV + 3y = , 
y'(o) = o, 

V" (oo) = 1 . 

The connection parameter associated to that problem is: 

F(0)=7, 

and one must determine the value 7 = 7* which corresponds to V* (x) 

^^Not always because the maximal value of _R may be larger than the radius of convergence of the series (see section 3.5.3.2). 
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(20) 
(21) 
(22) 



(23) 



It is convenient to work with the derivative V' (x) and to perform the foUowing ehange of variable: 



fix) = r{x) , 

f {x)=xg {x^) , 

then the system (20-22) is equivalent to: 

4:zg" -Azgg' -2g^ + {6- z) .g' + 2.g = 0, (24) 

g(+oo) = l, (25) 

with the new connection parameter defined as 

fc = 5(0) = -37. (26) 

Notice that, compared to the original BVP (20-22), there is no explicit initial condition (e.g. on g' (0)), 
it is replaced by the following condition of consistency on (24) at z = [a regular singular point of (24)]: 

which is automatically fulfilled when g (z) is expressed as a Maclaurin series the first terms of which are: 

giz)^k + i(fc -l)kz + ^{k - l)fc(8fc - 1) + ^(fc - l)fc2(34fc - 13) + O (z*) . 

It is relatively easy to verify the following (In appendix A we show how to obtain expansions like (27, 
28).: 

(i) eq. (24) admits two trivial solutions analytic in 2?+: 

(a) gl (z) = 

(b) 9*2 (z) = 1 

(ii) eq. (24) admits, locally to zq, the following singular solution (for a discussion of the reasons which lead 
to consider such expansions as solution of an ODE, see, e.g. [26,27]): 

5si„g (z) = + + - + ^ [l°g - ^ol) ' ^ ^o)]> ' (27) 

(Z - Zq) 4zo 

which depends on the two arbitrary constants zq and C (movable singularity) . Provided the expansion 
(27) converges, it represents locally the general solution of (24) which, thus, has necessarily a domain 
of definition smaller than D"*" except, perhaps, when zq — > +oo. Consequently, g* (z) does not belong 
to that two-parameter family, it is a particular solution of (24). 

(iii) eq. (24) admits a particular (one-parameter) solution which, asymptotically for large z, behaves as: 

g{z) ^ l + G^-2/5 + ^^-4/5 + o^-6/5\ (28) 

in which G is an arbitrary constant. 

The asymptotic one-parameter solution (28) satisfying (25), indicates the possible existence of the solution 
g* (z) as the envelope of the two-parameter family of locally singular solutions (27). By reversing the BVP 
(starting from the infinite boundary and with the connection parameter G), it is likely that a solution 
associated to any value of G will display a singularity somewhere before reaching the origin, except for a 
particular value G* (the envelope) which corresponds to g* (z). 

In conclusion, if we except the two trivial solutions g\ (z) and g2 {z) , the solution g* (z) is the only solution 
defined in In practice it is the envelope (zq — ?> -f-oo) of the two-parameter family of locally singular 
solutions (27) with C (continuously) adapted to the initial conditions. 

Since it is not singular on the positive real axis, one may expect that the domain of analyticity of g* (z) 
is sufficiently large to allow an analytic continuation of its Maclaurin series from the origin up to infinity. 
Moreover, since it belongs to a discrete set of nonsingular solutions [i.e., the set {gl (z),(72 {z)Tg* [z)}], 
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"minimal" procedures should work. This is confirmed a posteriori with the help of Padc approximants and 
conformal mappings. 

With 50 terms in the Maclaurin series (M = 50) k* has been estimated by Amore and Fernandez in [11] 
using the Pade-Hankel method yielding 20 digits of accuracy : 

k* w -0.22859820243702192438 . (29) 

Since, in a "minimal" procedure like the Pade-Hankel method, the condition at infinity is not imposed 
explicitly, it is almost indifferent whether one considers a basis of diagonal or a basis of sub diagonal 
approximants, in addition the two expected trivial solutions are also determined by observing (among the 
numerous values of k proposed as M grows) the presence of two additional convergent series of zeros towards 
and 1. 

According to the considerations of section 2.3.1.3 the success of the Pade-Hankel method is most certainly 
related to the possibility of moving to infinity (zq — > oo) of the movable singularity (27). But the broad 
accuracy of (29) is also surely a consequence of some particular analyticity properties of g* (z). Let us look 
at them. 

Knowing the value (29), even with fewer digits, one may explicitly use Pade approximants (diagonal or 
sub diagonal) on the Maclaurin series of g* (z) (with k ~ k*) with a view to determine the singularity 
distribution of g* (z). The procedure is not new (e.g., see Andersen and Geer [25]) and may be described as 
follows. The denominator of the rational function constructed from the series has some complex zeros which 
often reproduce the singular points of the actual function that it approximates. In general such singular 
points appear as the zeros (of the denominator) which remain stable as the order of the series grows (the 
spurious zeros being unstable). When the function to be approximated has a cut or a branch point it is 
simulated by an accumulation of poles on the cut. This is precisely what happens in the present case. For 
k sa —0.228598, the zeros of the denominator of the Pade approximants of a Maclaurin series (truncated to 
some varying value of M) accumulate exclusively on the negative side of the real axis of z. The stable zero 
(under changing the order M of the series) the closest to the origin controls the radius of convergence and 
gives R — 5.7217. Consequently, the angular sector of apparent analyticity would be the whole plane cut on 
the negative real axis and this corresponds to setting a = 2 in eq. (15). The values: 

i? = 5.7217, a = 2, (30) 

are thus the theoretically optimal values for the two parameters of the conformal mapping method. The 
latter value of R corroborates that obtained previously by using the d'Alembert rule [3] on the terms of the 
Maclaurin series for k k, k* . 

In ref. [24], the conformal mapping method has provided a better estimate of k* than (29) using values 
of R and a close to the maximal values (30) and with M = 120: 

k* w -0.228598202437021924373656107397047188424054430086. (31) 

It is important to realize that, though one does not know a priori the exact value of fc*, this (high) accuracy 
is internally controlled by the method itself. Indeed the convergence of the method may only be spoiled by 
ill controlled singularities. If the convergence is manifest, then the validity of the assumed properties of 
analyticity is confirmed and in turn the number of stable digits is secured. 

With the present case of the Polchinski ODE (24), the Pade-Hankel method is slightly more efficient 
than the conformal mapping method for relatively low orders. However it is very time consuming due to 
the necessity of calculating larger and larger determinants as the order M grows. The conformal mapping 
method is a bit less rapidly convergent but finally takes over because the calculations may be easily pushed 
to higher orders. 

Let us mention also that with the conformal mapping method, both procedures ("explicit" and "minimal" , 
see section 2.3) yield roughly the same accuracy. As for the Pade method, though the behavior (28) at the 

^*The authors had estimated k* — 1 instead of fc*. 

When the order of the series is not very high, the maximal values of R and a are not necessarily the values which give the 
best apparent convergence of the conformal mapping method. One observes, however, a tendency to favor the maximal values 
at very high order. 
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second boundary of the solution cannot be reproduced by a rational function one can impose an explicit 
condition. Actually, due to the condition (25), one can obtain an estimate of k* by using diagonal Pade 
approximants to represent the solution and by merely imposing that the ratio of the highest degrees of 
the respective polynomials P and Q be equal to 1. But considering the great efficiency of the "minimal" 
procedure, the recourse to the "explicit" version of the Pade method is not pertinent here. 

Let us emphasize that once k* has been determined within some accuracy, each of the two methods (except 
sometimes accidentally with a Pade approximant) provides an explicit analytical (approximate) expression 
for the solution: a rational fraction (a Pade approximant) or a convergent series in powers of w (z) as given 
by (15). This remark is obviously valid in any occasion where the methods work. 



3.2. The Thomas-Fermi equation for the neutral atom 

The Thomas-Fermi problem for the neutral atom [17,18] writes in terms of a function u (x): 

u" - \/— = 0, 
V X 

ii(0) = l, 
u (+oo) = , 

with the connection parameter k = u' (0). 

As in the previous example, it is convenient to perform a change of function and of variable: 



g{z) = Vm^' 

the original problem is then equivalent to: 

z[gg"+g'^]-9g'-2z^g^^0, (32) 
g (0) = 1 , (33) 
g{+^)^0, (34) 

and the connection parameter is: 

k = g" (0) . 

The first terms of the Maclaurin series of g (z) are: 

, , kz^ 2z^ Pz'^ 2kz^ 1 f^k^ \ efc^z^ ^, o, 

^(^)-^ + - + --^-^+18(--V'^^''''^'^- 
The Thomas-Fermi equation has much more been mathematically studied than the previous one. In 
particular it is known, by a theorem of Mambriani [28] (see [29]) that g* (z) is unique. Sommerfeld [30] has 
calculated its behavior for large z, and one may easily verify (see appendix A) that (32) admits the following 
asymptotic expansion as a one-parameter solution: 

g{z) ^ + £ + s = x/73-4, (35) 

z-i.-1-oo z'^ Z^ 

in which G is the only arbitrary constant. As in the Polchinski case, it exists a value G* which should 
correspond to g* (z) and, presumably, it is the envelope of a family of singular solutions satisfying the initial 
condition (33). 

Indeed one may easily show that, similarly to (27) for the Polchinski case (see appendix A), the general 
solution of (32) has a two-parameter movable singularities which, locally to any zq, behaves as: 
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5 1 20 1 35 1015 , 

1175 , ,o 43805 , ,o 

104976zo5 ^ 157464zo^ 
136949335 4 2338092205 5 

"9'^«n«c^.^R«.„7(^~^0^ + 91/1977011 9. „8 ~ 



r Z ~ Zn H 7r[Z — Zn) 

104976zo5 ^ 157464zo6 ^ 

136949335 4 2338092205 

"238085568V ''' 2142770112^^^" 

189476597645 g 37136484275 

"77139724032V ° ''" 4132485216V" 

+(z - zo)« f C + 

^ ^' \ /l9fi1fi9f^'?7Q^ll / 



- Zo) 



i7 



MOU»»UU10gUZo-Z|J \ 

4261625379zii / " ' 

in which zq and C are two arbitrary constants . As in the previous case, it represents locahy the general 
solution of (32) which, thus, has necessarily a domain of definition smaller than T)'^ except, perhaps, when 
zo — ?> 00. 

Recalling section 2.3.1.3, we are not surprised that the Pade-Hankel method provides an easy determi- 
nation of k* . With A/ « 60 Amore and Fernandez [11] (see also Fernandez [31,32]) using this method have 
obtained the following estimate: 

k* « -1.58807102261137532 . (37) 

In the circumstances and as in the Polchinski case of section 3.1, with regards to the efficiency of the 
"minimal" procedure, it is not necessary to make use of the "explicit" procedure with Pade approximants 
to determine k* . But with the knowledge of fc*, one may construct a rational fraction (only once) to obtain 
an explicit (approximate) solution. As explained in section 3.1, that approximate explicit solution may be 
used to roughly determine the distribution of the singularities of g* (z). 

For k K, —1.5880710 and values of M up to 120, we observe, compared to the Polchinski case, a wider 
dispersions of the stable singularities, but they remain confined to Re (z) < with a cut on the real axis 
starting at z = —0.6658. The maximal angular sector of analyticity appears to be characterized by the value 
a = 1.1. Hence the values: 

i? = 0.6658, a = 1.10, (38) 

arc the theoretically maximal values for the two parameters of the conformal mapping method. This corre- 
sponds to a smaller range of analyticity than in the Polchinski case. 

According to the remark of footnote 15, for practical efficiency, we have chosen R = 0.65 and the convenient 
a = 1 which are close to the maximal values (38) and with M = 200 we have obtained the following estimate 
of k*: 

k* « -1.5880710226113753127189 ± 7 x 10"^^ , 

which is only slightly better than the result (37) previously obtained by using the Pade-Hankel method. 
Actually, at low orders, the calculations clearly show a greater efficiency of the Pade-Hankel method com- 
pared to the relatively poor results of the conformal mapping method (which provides a better estimate 
because high orders may be treated). 

We observe also a smaller rate of convergence of the two methods, compared to the Polchinski case. The 
effective value of a being smaller than in the this latter case (where the maximal allowed values a = 2 is 
reached), one may conclude that the (relative) greater difficulties encountered with the Thomas- Fermi ODE, 
is a direct consequence of a less favorable distribution of singularities. 



^®The late presence of the logarithm in this expansion prevents the ODE from having the Painleve property (for a review see, 
e.g. [26,27]). Notice also that the explicit calculation of all the terms of (36) to reveal the presence of the arbitrary constant C 
in front of the monomial (z — 20)*, is not necessary. Such a term is sometimes named a resonance. Once the leading term has 
been determined, it is sufficient to find the value of the degree n of an arbitrary additional monomial which linearly contributes 
to zero in the ODE. In the present case, one finds n = 8. (See appendix A.) 

^7 Compared to the Polchinski case for which (31) is much more accurate than (29) whereas M = 120. 
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3.3. The third grade fluid in a porous half space 



In this section we consider the equation of the flow of a third grade fluid in a porous half space which has 
been investigated by Hayat et al [19] and Ahmad [20]. The formulation of this problem in terms of / {z) is 
as: 



/" + &i/'Y'-^/r-c/=0, (39) 

/(O) = 1, (40) 
/(+^) = 0. (41) 



in which 6i and c are two constants. 
Using a rescaling of the form: 



lh\c ' 



the problem (39-41) may be written as follows: 

1 
3 



5" + 5'V'-^.95"-,9 = 0, (42) 



g{0) = A = ^c, (43) 
g{+^) = 0. (44) 

The connection parameter is again: 

k = g'{0) . 

For information, the particular values bi = 0.6 and c = 0.5 considered in [20] had yielded the following 
estimate (using a shooting method [20]): 

/*' (0) w -0.678301 , 

which, in our conventions, corresponds to 

A = ^/hc ~ 0.5477225575 , (45) 
k*^\/fhf*' (0) w -0.525410. 

Having fixed the values of bi and c to the values indicated above, the first terms of the Maclaurin series of 
g {z) read as follows: 

T (k^ + 3) z2 k (fc2 + 3) l5k^ + lOk^ + 3) 

\-kz + — ^= \ !^ ^-^ = ^ h O (zn . 

10 2V30(fc2 + l) 90(fc2 + l)^ 

It is easy to verify that the ODE (42) does not admit, locally to zq, any solution of the form^^ : 

ffsing (z) ^ A {z- zqY , (46) 

z— >-zo 

but for s ~ 1 which makes gsing (z) analytic. 

Actually the general solution has a movable essential singularity at infinity since for large z, (42) asymp- 
totically admits the following two-parameter family of behavior as solution (sec appendix A): 

g{z) ~ Cie75 - ( ^z + C2 I e^Tf -^0(ze"^7f), (47) 

in which Ci and C2 are two arbitrary constants. Notice that the limit Ci — )• in (47) is singular, this 
indicates the possible existence of another kind of solution. 



Since the ODE is autonomous, it is sufficient to look for a solution of the form A. z" (because of the translation invariance 
under the change z — > z + zi). 
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Actually, there is a one-parameter family of solutions of (42) which, for large z > 0, behaves as: 

giz) ^ Ge-_^e-3^ (48) 

This kind of solution (which goes to zero at the second boundary) should correspond to g* (z) for a 
particular value G* . It is the envelope of the solutions which, for z +oo go to ±oo according to the sign 
of Ci in (47). Indeed it corresponds to the singular limit Ci — > of the general solution corresponding to 
(47). 

Despite the probable uniqueness of g* (z) as a solution of (42) defined in 2?+, the Pade-Hankel method 
works with difficulty yielding the following poor estimate of k*: 

A:* w -0.525, (49) 

with the order of the Maclaurin series limited (to save time) to M = 17. This poor result is in agreement 
with the absence of movable pole and confirms the considerations of section 2.3.1.3 on the nature of the 
Pade-Hankel method. 

The "explicit" procedure with diagonal Fade approximants has not improved the result. 

One may associate this poor behavior with the fact that essential singularities located at infinity seem to 
be analytic when they are seen from the origin, they are more difficult to detect than poles or cuts located 
at a finite distance from the origin. 

Despite the poor accuracy of (49), Fade approximants may be used to roughly give the distribution of 
singularities of g* (z) in the complex plane of z (see figure 1). 

This is sufficient to determine values for the two parameters of the conformal mapping method: 

1.59810, a = 0.547. (50) 

However, because the first real singularity is far from the origin (see figure 1), the largest possible value of 
the parameter R of the mapping method does not coincide with the radius Rq = 1.59810 of convergence of 
the Maclaurin series. This is true each time the closest singularity is not located on the negative real axis (the 
maximal value of R is always associated with the closest negative real singularity). But the largest possible 
value of R is not necessarily its optimal value as parameter in the conformal mapping method because it 
induces a decrease of a compare to the choice R = Rq. This is why we have kept the values (50) with which 
we have obtained the following estimate of k* (for M = 80): 

k* = -0.52541018 ± 1 X IQ-^ . 

We have observed that the "explicit" mapping procedure has appeared more efficient than the "minimal" 
one which has given a slightly less accurate estimate of k*. This effect may be associated with the remark 
done above concerning the greater difficulty of detecting an essential singularity located at infinity compared 
to a pole or a cut. We note also that despite a less favorable distribution of singularities the mapping method 
has, this time, clearly given a better estimate than the estimate (49) obtained with the Fadc-Hankel method 
(confirming again the nature of this method, see section 2.3.1.3). 

Actually, there is a much more efficient method to solve our problem. It is based on the following change 
of variable as proposed in [20]: 

w = , 

which transforms the original problem (42-44) into: 



g"{w) [l + w'^ g' [wf] +wg'{w) [l + w^g'{w 
-g{w 



l + lw'g'iw)' 




Actually there is also the same kind of solution defined in 'D~ that corresponds to the value k = —k* and its presence is well 
observed as a convergent series of zeros in the Pade-Hankel method. This solution is not observed with the "minimal" version 
of the mapping method since only 1)+ is actually mapped onto the unit circle. 
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Fig. 1. Singularities (small circles) in the complex plane of z of the solution of the ODE (42) for k ft: —0.52541 and M = 51. 
The three solid small circles indicate the position of the most important stable singularities of the Fade approximants. The two 
open circles on the real axis are spurious (unstable) singularities. The fact that the stable (black circle) real singularity is far 
from the origin allows a freedom in the choice of the parameters R and a of the conformal mapping method and in particular 
R may take on a value larger than the radius Rq of convergence of the Maclaurin series, see text for more details. 



with a new connection parameter: 

g' {Q)^G. 

On expressing the generic solution g (w) = g[z {w)] as a Maclaurin series about w = 0, the first terms of 
this expansion read [see (48)]: 

giw)=Gw- -G' + —G' - G^w' + G^w' + 0( . (51) 

' 12 432 15552 103680 ^ ' ^ ' 

It is a matter of fact that the (simple) sum of the series (51) converges in the range w g [0, 1]. Then by 
imposing that this simple sum satisfies the required condition at the second boundary {w = 1), namely: 
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r (1) = g* [0] = ^Jbic = 

we have very easily obtained, with 60 terms in the series, the following estimate: 

G* = 0.56044886606934678 ± 1 x 10"" . (52) 

Using the relation: 

g' (z) = -wg' [w) , 

it comes: 

k^--g'{l) = -G+\G^-^G^ + --- , 

which, for G = G* as expressed in (52), finally gives the following estimate for the quantity originally of 
interest: 

k* w -0.525410175091336. 

Using Padc approximants on the Maclaurin series in powers of w (51), we have verified that, for G ~ G* 
as given by (52), the singularities of g* (w) are disposed in the complex plane of the variable w in such a 
manner that the radius of convergence of the series is actually larger than one. Consequently, for the value 
of A considered [sec cqs. (43, 45)], the series actually converges and the recourse to an analytic continuation, 
such as a Pade approximant, as done in [20], is not necessary to obtain an explicit analytical form of the 
solution. However, when A grows, this radius of convergence decreases and becomes smaller than one for 
some ^ < 1 so that the simple sum no longer converges at w = 1. It is very likely that for those values of A, 
the use of one of the two methods considered in this paper and applied to the series in powers of w, would 
have some efficiency. 




3.4. The Falkner-Skan flow equation 



In this section we consider the boundary layer Falkner-Skan equation for wedge. By choosing an especial 
case of magnetic field, the boundary layer similarity equation obtained is (Abbasbandy and Hayat [33]): 

g'" + gg" + (i - .9") - (g' - 1) = , (53) 

5(0) = 0, ff'(0) = 0, (54) 

g' (00) = 1 . (55) 

The connection parameter is defined as: 

k = g" (0) . 

The Falkner-Skan ODE (53) is much more complicated than the previous examples. Its complete discussion 
is out of the scope of this article and, for the sake of our illustration, we limit ourselves to the case /3 = 4/3 
and A = 2 already studied in [33] and for which the first terms of the Maclaurin series of g [z) are: 

(z) = z + -z^ ^ [bk^ - 64) z'^ z^ H z^ + O (z^) . 

^ ^ ^ 2 9 6 360 ^ ' 135 2835 ^ ' 

One may verify that (53) admits the following three-parameter movable singularity as local solution for z 

close to zq (see appendix A): 

ffsing (z) ~ - r, h Gl (z - Zo)"+ + G2 (^ - ^o)"" , (56) 

p — Z Z — Zq 

_ 5/3 - 4 ± \/25/3^ - 16/?- 32 
2(/3-2) ' 

in which zq, Gi and G2 are arbitrary constants. 

For the expansion (56) to make sense, one must have Re {a±) > — 1, which occurs for: 
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/3 < f 2 - sVe) w -0.855755 or /S > 2 . 
25 V / 

These constraints exclude the value of interest /? = 4/3, for which we have: 

a± = -2±iV5, 

and, in the circumstances, the expansion (56) does not locally represent the general solution. 

There is, however, another movable singular solution which is of interest to us. It involves only two 
arbitrary parameters (thus it is not a general solution) , and locally about zq behaves as : 

{z - zof (140 - 27CiA') / 5 



+ 0(iz-z.r). (57) 

We observe that this expression is singular when Ci — > and as in section 3.3 [eq. (47)] this presumably 
indicates the existence of a unique solution defined in 2?+ (the envelope of the singular solutions) which 
could be g* (z). 

Eq.(53) admits also two particular exact solutions of no interest to us since they do not correspond to the 
required boundary conditions: 

gi {z) = z + Bi, 

92 {z) = Z + B2. 

These two linear behaviors appear also as leading terms of two asymptotic expansions (z — > 00). In fact, 
for the value of interest /3 — 4/3, (53) admits, asymptotically when z — >■ ±cx), a one-parameter solution 
which is of interest to us: 



g{z) ~ z + Gz-^~^^-^e- 

z— >-iboo 

1 



1 (aV2/3 + 2)(aV2/? + 3) ^ ^ 
2z2 ^ 



^G2(/3 - l)z-2^ -4^-5e-^ [1 + (z-2)] + o(^«e-3^) , (58) 

in which a has to be defined in terms of [3 and A. It is potentially the asymptotic behavior of the solution 
g* (z) for [see eq. (55)]. 

For some values of /3 and A (not for the value /? = 4/3 however) it exists another particular asymptotic 
solution (corresponding to the second exact linear solution) that behaves as: 

g{z) — z, 

with exponential corrections. This behavior does not satisfy the condition (55) at the second boundary and 
does not exist for /3 = 4/3. We shall not discuss it further. 

The existence of the one-parameter asymptotic solution satisfying the required condition (55) at the 
second boundary indicates that g* (z). if it exists, is likely to be unique. It is probably the envelope of the 
family of the two-parameter solutions having locally the movable singularity (57). 

It is thus not surprising that the "minimal" Pade-Hankel method succeeds in determining a unique value 
of k* [33] (with M ~ 30). We have redone the calculation for a greater M = 51 and obtain the following 
improved estimate (for other values of the parameters /3 and A, see [33]): 

k* = 2.43949894 ± 1 x 10"^ . (59) 



Since this is a local expression of a particular solution, it may well have only poles (no logarithm) without implying that the 
Painleve property is satisfied. 

^'^For /3 < 0, there is an additional arbitrary correction term to (58) of the form Gz^P~^^. 
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Fig. 2. Singularities (small circles) in the complex plane of z of the solution of the Falkner-Skan ODE (53) with A = 2 and 
/3 = 4/3, for fc 2.43949 and M = 100. The four solid small circles indicate the position of the most important stable 
singularities of the Fade approximants. The values i? = 2.28 and a = 0.51 arc the theoretically optimal values of the parameter 
of the conformal mapping eq. (15), see text for more details. 

Using again Pade approximants we may have a look at the distribution of the singularities of g* {z) in the 
complex plane of z, see figure 2. From this distribution we obtain: 

i?w2.28, aw 0.51. 

Using values of i? and a close to these estimates, the conformal mapping method yields convergent 
procedures. However, even with M ~ 130, the result does not actually improve the estimate (59) obtained 
by the Pade-Hankel method with M = 51 (see fig. 3). The better efficiency of this latter method compared 
to the primer one is reestablished and reinforced due to the hardening of the analyticity properties of the 
solution. 
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Fig. 3. [Falkner-Skan ODE (53) with A = 2 and /3 = 4/3] Comparison of the evolutions with M of the estimate of k* according 
to the method used. Two configurations of the conformal mapping for a fixed a = 1/2 (resp. _R = 2.3 and -R = 2) are shown 
(see text for more details). 

3.5. The Blasius Equation 



3.5.1. Generalities 

The Blasius equation is the Falkner-Skan equation (53) with /? = A = [1]: 

.9"' + m"-0. 



(60) 
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Let us first consider the general BVP of finding a solution to (60) which satisfies the following conditions: 



.9(0) = a, g'{0)^b, (61) 
g'(+c^)=i?, (62) 

and the connection parameter k is defined by: 

k = g" (0) . (63) 

Eq. (60) may be integrated to give: 

g" (z) ^ k cxp £ g{t)^ , (64) 

consequently g (z) is either convex if fc > 0, or concave if fc < or linear if fc = 0. 

For the sake of our illustration, we limit ourselves to the case B = 1 with b E [0,1) which corresponds to 
convex solutions (for a discussion of the concave solutions, see for example [34]). 

There are two explicit exact particular solutions to (60), the first one corresponds to the linear case: 

g\in (z) ^ a + bz, 
and the second one has a movable singularity: 

3 



(65) 



Z~ Zq 

corresponding to the particular conditions at the origin: 

- 3 ,3 6 

Oi - — - , - — ■ — 2 7 ~" — ^ • 

Zq Zq 

Because the ODE is autonomous, these two kinds of particular solutions (linear and singular) appear 
again, either as asymptotic limits when z — > +oo, or as singular limits in the two main kinds of solution 
(concave and convex): 

- Similarly to the Falkner-Skan case, there is a three-parameter movable singularity representing locally the 
general solution [see eq.(56) for /3 = 0] 

ffsing [z) ^ 3— h Ci (z - zo)"^ + C2 (z - zo)"~ , (66) 

Z Zq 

a± = l± iV2, 

which occurs only in the concave case. 

- For large z, two asymptotic solutions exist: 

(i) a one- parameter asymptotic solution (Gi is arbitrary): 

ffasy (z) ^ - + % + ^+ (z-4) , 
z->oo z Z^ 3Z'' ^ ' 

which does not correspond to the convex case. 

(ii) a two-parameter asymptotic solution (Gi and G2 are arbitrary): 



G2 

ge^syiz) ~ Giz+^e —, (67) 

which is compatible with the convex case if Gi > and candidate to be the solution of interest if Gi = B = 1. 

Notice that (67) replaces the one-parameter asymptotic solution (58) of the Falkner-Skan case. This is in 
agreement with the disappearance of the particular two-parameter local solution with movable singularity 
(57) and, consequently, of their envelope and the possibility of having an isolated (unique) solution defined 
in 2?+ (at least in the convex case). Let us confirm this by the following considerations. 
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In the convex case that we consider, we have g (t) > a + bt and from (64) we can write: 

< g"iz) < ke-"'-^^'^ , 

consequently the domain of definition of these solutions is T)^ and there is a continuum of such solutions 
corresponding to any /c > 0. Of course only one of them gives g' (+oo) = 1, but if that latter condition is 
not explicitly imposed, it is unlikely that the right value of k* be determined. Consequently, the "minimal" 
procedures are inefficient here. This is the reason why the Pade-Hankel method has failed to solve the Blasius 
problem [11]. 



3.5.2. The Blasius problem 



The Blasius problem consists in finding the solution of (60) which satisfies the following boundary condi- 
tions (i? = 1, a = 6 = 0) [let us mention the interesting Boyd's review on this problem [35]]: 

.9(0)=0, .g'(0) = 0, (68) 
g' (c«) = 1 . (69) 

With the conditions (68) at the origin, the function ,9(2) appears to be of the form: 

g (z) = z^u (z3) . (70) 

Indeed Blasius [1] has given the general expression of the Maclaurin series of g (z) corresponding to the 
initial conditions (68), it reads: 

(-1)^. 

{1 i = Qui = 1 
^ (3^-1)! . . o ' (^2) 

^ (3j)! (3* - 3j - 1)! ^ ' 

which corresponds to the following first explicit terms: 

„ {-A = -z^ - —z^ + ii^z« ^^z" + ^^^^^^ z" + O fz") 

^ ^ ^ 2 120 40320 532224 29059430400 ^ ^ ' 

Moreover the particular form of (60)] induces the following scale invariance [36]: 

(73) 

5-> 5 = cr5(CTz) , (74) 



z 

z — > z = — 
a 



so that, without changing the conditions (68), the new connection parameter k = g" (0) = <j^k may be set 
a priori equal to any value, say qo, in the Maclaurin series of the Blasius problem, i.e.: 

k^q,=^a=[^) , 
in which case, the value of g' (00) at infinity reads: 



ff'(oo)=(|j ' <?'(«.). 

Hence, for the required condition (69) at infinity, we simply read the value of k* as the result of a simple 
Cauchy problem with the conditions (68) and g" (0) = q^: 

k* =q^[g'{^)]-^'\ (75) 



2/3 
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Wc have thus integrated numericahy (60) for fc = 1 (using Mathematica) and obtained (see also Boyd 

[13]) 

k* w 0.469599988361013304509 . (76) 

We may also sum the Maclaurin series using a Pade approximant or an adequate conformal mapping. To 
this end, it is useful to perform the change of function 5 — >■ u of (70) and to work with the new independent 
variable 

X = . 

If we note u (x) the function corresponding to g (z) introduced above with g" (0) = fc = 1, then we have: 

[g' {z)f = X [2u (x) + 3xu' {x)f . (77) 

One may then perform a diagonal Pade approximant of the corresponding truncated series and obtain an 
estimate of k* by considering the following limit of the rational fraction: 

fc* = I lim Pade (a; [2u {x) + 3xu' {x)f^ | . (78) 

This way we have obtained (with AI = 200): 

k* « 0.46959998836100 , (79) 

but the observed convergence is not limpid. 

We may also extract, from the Pade approximant of u (x) with a given fc, the corresponding distribution 
of singularities in the complex plane of x. For the value k = 0.4696 close to (79), we obtain a distribution 
of singularities located in the Ihs half plane with a cut on the real axis starting at a; = —65.133 and some 
singularities elsewhere which induces the following values for the parameters of the conformal mapping: 

i? = 65.133, a = 1.57. 

As noted in footnote 22, this value of R agrees with the preceding estimate of [13] in which reference 
one may find a discussion of the distribution of the singularities in a particular choice of normalization for 
which a = 3/2. For other normalizations and other values of k, the estimates for a and R change. We have 
performed a conformal mapping on a series with M ~ 590 for a fixed fc = 1 and obtained after a single sum 
of the scries [using (75) with qq ~ 1]: 

k* w 0.46959998836102 . 

That estimate is not more accurate than (79), showing again the great efficiency of the Pade approximants 
(though it would be difficult to improve the result whereas this would be relatively easy with conformal 
mappings) . 

For illustrative purposes we present Table 1 in which we present our results compared to various estimates 
of k* encountered in the literature. 

The Blasius problem is finally easy to solve because of the conjunction of the scaling invariance and of 
the initial conditions (68) which induce the effective independent variable x = z^ and increases by a factor 
three the angle of the domain of analyticity of the solution. This favorable conjonction no longer exists with 
different initial conditions as shown in the following subsection. 



There are several normalizations used in the literature. The most common is to choose a function / (z) = 2g (2) with B = 1 
and the connection parameter kg = f" (0), in which case one has = 0.33205733621519630 [13]. For the choice g (2) we have 
made, we get k* = k^ if one fixes B = 2 (the genuine Blasius problem). With the present choice of -B = 1, using the scale 
invariant properties (74), we have k* = ^/2kg which corresponds to (76). As for the radius of convergence Ro = 5.6900380545 
of the Maclaurin series in powers of 2 given in [13], it is related to our R = 65.113 for the series in powers of x = 2"^ via the 

relation R = (_Ro/V2)^ = 65.13291560. 
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Source 


fc* 


Mcth 


present work 


0.469599988361013304509 Integ 


Boyd [13] 






present work 


0.46959998836102 


MAF (590j 


present work 


0.46959998836100 


PAD (200) 


Fazio [37] 


0.469599988361 


SHO 


Asaithambi [38 


0.46959998847 


SHO 


Parlange [39] 


0.46959999 


Integ 


Yu-Chen [40] 


0.4695997 


DTM 


Liao [41] 


0.469500 


HAM 


Howarth [42] 


0.46960 


SHO 


Blasius [1] 


0.4695 


Matching 



Table 1 

(Blasius problem) Estimates of k* of the present study compared to other estimates encountered in the literature. The last 
column on the right indicates the method used. "Matching" : matching of the sum of the Maclaurin series with the sum of an 
asymptotic expansion; "SHO": shooting method (Cauchy+Ncwton); "Integ": a single integration (Cauchy); "HAM": homotopy 
analysis method; "DTM": Differential Transformation Method; "PAD (200)": Pad approximants on the Maclaurin series with 
M = 200 using eq.(78); "MAP (590)": a conformal mapping on the Maclaurin series with M = 590 using eq.(75). 

3.5.3. A more general convex solution 

Let us consider the BVP of finding a solution to (60) which satisfies the foUowing conditions (a similar 
problem has been studied by Allan in [43]): 

g(0) = 0, 5'(0)=6, (80) 
5'(+cx3) = l, (81) 

as explained in section 3.5.1, the corresponding solution is convex provided that 6 < 1, in which case the 
associated connection parameter k defined in (63) takes on a positive (still unknown) value /cq . 



3.5.3.1. A simple dichotomy algorithm For the sake of our illustration, we first consider the case b = 1/4. 

Before considering any series expansion, let us show an interesting aspect of the invariance of (60) under 
the rescaling (74, 73). Though the BVP could no longer be transformed into a simple Cauchy problem, we 
may manage an easy dichotomy procedure to determine feg. 

Denoting by [cr = 1 : {b, k^] 1}] the original BVP corresponding to the boundary conditions (80, 81), then 
by virtue of the scaling invariance, the following equivalence stands: 



1 : {&, fc*; 1}] ^[a= {k*)'/' : {b (k*)-'^' , 1; (k*)-'^'] 



Consequently, to determine the value fcg (associated to b), wc may proceed as follows. Choosing an arbitrary 
number bo, a single integration of the ODE (60) with the initial condition g (0) = 0, g' (0) ~ bo, g" (0) ~ 1 
provides the value g' (oo) = Bq which should be equal to b^/b otherwise one tries the new value g' (0) = 
bo+bBo gQ g^j^jj gQ fQj-|;}j until g' {oo) « g' (0) /b, which finally gives k^ « [g' (oo)]"'^^^. Of course the 

speed of the method depends on whether the initial value bo is close or not to b (fcg)"^^'^. By this dichotomy 
procedure we have obtained the following estimate for b ~ 1/4: 

k* = 0.429540737652735 ± 1 x 10"^^ . 

At this stage, we already mention that the Taylor series methods, which we are presently interested in, have 
failed in providing any reliable estimate of fcg while a [10,10]-Homotopy Pade method (for an introduction 
see [44]) gives, with a relative easiness, the estimate fcg = 0.429543. Let us look at the reasons why the 
Taylor-series-based methods are not efficient there. 
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3.5.3.2. Taylor series and the analyticity properties of a solution Wc already know that the "minimal" 
procedures (in particular the Padc-Hankel method) cannot work because of a continuum of solutions defined 
in 2?+ (sec section 3.5.1). In addition the singularities in the complex plane of the independent variable z 
of the solution cannot be moved toward the left-hand side as in the Blasius problem. Consequently some of 
them are located in the region Re (z) > in such a way to limit the efficiency of the Pade and conformal 
mapping methods. With a view to an illustration of the origin and the extent of the difficulties encountered, 
we give up solving the EVP and turn our attention to the Cauchy problem associated to the following initial 
values: 

.9(0) =0, g'(0) = l/2, g"(0) = l, 
for which a single integration provides the asymptotic value: 

g' (oo) w 1.78283921978496662798847022654 . (82) 
The first terms of the corresponding Maclaurin series read: 

g{z)^- + -^-- — + — + ^^z' + ^^z»-^^z' + 0(z'°) . (83) 
■^^ ^ 2 2 48 120 960 10080 322560 483840 ^ ^ ^ ' 

By performing diagonal Pade approximants on this series, we obtain, the distribution of singularities as 

shown in fig. 4. From this we observe that: 

(i) the largest possible value of the parameter R of the mapping method does not coincide with the radius 
of convergence of the Maclaurin series (noted Rq in fig. 4). But, this time the radius of convergence Ro 
cannot be chosen as reference to determine the optimal value of a as it would follow from the idealized 
scheme of fig. 1 of [24] , this is because: 

(ii) farther singularities in the half plane Re (z) > would prevent the convergence of the mapping method, 
so that 

(iii) the effective maximal optimal value of a (0.31) is small and can be significantly enlarged only if R is 
considerably diminished. 

As consequence of this unfavorable singularity distribution, the mapping method for the optimal values 
R = 3.8 and a = 0.3, converges only very slowly (sec fig. 5) and, with M ~ 350, we have obtained the 
poor estimate g' (oo) = 1.78 ± 0.06. Smaller values of R yield slightly improved convergences which remain 
slow however. Nevertheless the convergence (even slow) of the method indicates that the distribution of 
singularities of fig. 4 is correct. Despite this fact, Pade approximants do not yield any estimate of g' {oo). 

4. Summary and conclusion 

Wc have presented, compared and discussed the efficiency of three quasi-analytic methods for solving a 
BVP: the Taylor series, the Pade and the conformal mapping methods. After having reminded that the first 
method almost always requires the recourse to analytical continuation procedures to be efficient, our interest 
has been focused on the latter two. 

We have emphasized that the Pade-Hankel method is indeed the "minimal" version of the standard Pade 
method and explained why it is not always efficient. In particular we have shown that its efficiency is closely 
related with the possibility of pushing to infinity a movable singularity. The conformal mapping always works 
but (obviously by construction) its efficiency depends on the analytic properties of the solutions sought. 

With a view to present a fairly comprehensive discussion, we have successively considered six different 
configurations of ODE and initial conditions which yield more and more constraining analytic properties. 
We have shown that Pade approximants may be used to determine correctly the singularities distribution 
which allows an efficient determination of the values of the two parameters R and a of the conformal 
mapping method [sec cq. (15)]. In each case of ODE considered wc have empirically discussed the conditions 
of uniqueness of the solution defined in 2?+ = [0,+oo), this unique solution appears as the envelope of 
a family of solutions having moving singularities, when such a family of solutions does not exist then an 
infinite number of solutions are defined in and "minimal" procedures (like the Pade-Hankel method) do 
not work. 
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Fig. 4. Singularities (small circles) in the complex plane of z of the solution of the Blasius ODE for the initial conditions (80) 
with h = 1/2 and g" (0) = 1 (from the Padc approximant [35,35], discarding only some spurious singularities). The small solid 
circles indicate the positions of the stable singularities the closest to the origin. The large circle of radius Rq = 2.8616 indicates 
the domain of convergence of the Maclaurin series (83) in powers of z. The theoretical optimal values of the parameters of the 
mapping method are rather R = 3.8309 and a = 0.31. Farther on the rhs, singularities (not yet completely stable) may perturb 
the convergence of the mapped series. 



Contrary to the mapping method, the Padc-Hankel method (when it works) is comparatively more and 
more efficient as the analytic properties harden but it is limited by its awkwardness. The mapping method 
is more sophisticated and smoother. Provided the analytic properties of the solutions are favorable, it is 
useful when a high accuracy is required. 



Acknowledgments We thank B. Boisseau and H. Giacomini for their useful remarks and comments 
relative to this work. We are also indebted to the referees for their useful suggestions. 
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Fig. 5. [Cauchy problem for the Blasius ODE (60) with initial conditions a = 0, 6 = 1/2 and fc = 1] Comparison of the evolutions 
with M of the estimate of g' (oo) by the mapping method for a = 0.3 and two values of _R. The solid horizontal line indicates 
the expected value (82). The convergence is slow but manifest (see text for more details). 

Appendix A. Expansions of solutions about particular points 



Asymptotic expansions such as (28, 35, etc..) and expansions about a movable singularity such as (27, 
36, etc..) may be determined following a similar procedure. The main difference with the expression of the 
generic solution as a Taylor series is that the terms are not necessarily integer powers of the independent 
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variables. For example, it may involve also powers of logarithms. Sometimes the leading term is not a power 
but an exponential or a logarithm. This prevents us from using a general systematic algorithm to get the 
terms of such expansions. 

As an illustration we present below the example of the asymptotic expansion of the solution of the 
Polchinski fixed point equation. The obtention of the local expansion about a movable singularity is very 
similar and will not be made explicit here. 

A.l. Asymptotic expansion of the solution of the Polchinski equation 

Let us consider the ODE (24). We want to show that it admits the one-parameter solution (28) as z — > oo. 
Although not obliged, it is convenient to perform the following change: 

g{z) = h{y = l/z) (A.l) 

then (24) reads: 

4y3/^"(y) + 2y^h\y) + 4yh{y)h'iy) + yh' {y) - 2h{vf + 2h{v) = (A.2) 
Assuming that y is small, we try to find a local solution to this ODE under the following form: 

h^M-XV'' (A.3) 

where x S'lid 7 are the unknowns. It is easy to see that (A.3) generates three different powers of y in the 
ODE: y'^ , y'^^' and y^'^^ . To get a solution (locally valid) one must first determine the leading power for 
small y and second impose that its coefficient vanishes. These two conditions should determine the values 
of X and 7. 

Since y is small, the smallest power prevails, then we distinguish three possibilities: 

(i) 7 < 0, the leading power is 27, but the vanishing of the global coefficient imposes 7 = 1/2 (whatever 
x) what is incompatible with the hypothesis. 

(ii) 7 > 0, the power 7 prevails, but the vanishing of the global coefficient imposes 7 = —2 (whatever x) 
what is incompatible with the hypothesis. 

(iii) 7 = 0, the two powers 7 and 27 are equal and the vanishing of the global coefficient occurs for x = 
or X = 1- 

The leading term of the asymptotic expansion is thus 1. To get the next term we try the following form: 

h^y{y) = l + Xy'' (A.4) 

in which the second term is smaller than the first one, i.e. 7 > 0, for consistency. Introduced into the ODE, 
this expansion about y = gives the same set of powers of y as previously. The leading power is thus 7. The 
vanishing of the global coefficient implies that 7 = 2/5 whatever x- That is an acceptable solution. 

The arbitrariness of the amplitude provides the first arbitrary constant that is noted G in (28). 

The local expansion of the solution about y = may be continued by considering the next correction: 

/iasy(y) = l + Gy'/'' + xy^ 

and so on and so forth. 

In the present case, G is the only arbitrary constant appearing in the asymptotic expansion. It is important 
to know the maximal number of arbitrary parameters in the solution (sometimes called "resonances" in the 
expansion about a movable point). However, it may occur that the arbitrary constants appear lately in the 
expansion [e.g. see (36)] which, in addition, may be very complicated. There is a convenient way to quickly 
know the values of 7 that are associated to the arbitrary constants in the expansion without having to 
calculate it explicitly. Once the leading term of the expansion is known, one considers the generic correction 
as in (A.4) then it suffices to look at the vanishing of the linear-in-x contribution to the ODE. The resulting 
constraint provides an equation for 7. In the present example there is a unique acceptable solution to this 
equation: 7 = 2/5, there is thus only one arbitrary constant (G). 
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